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We study transport through a triple quantum dot in a triangular geometry with applied bias such 
that both singly- and doubly- charged states participate. We describe the formation of electronic 
dark states — coherent superpositions that block current flow — in the system, and focus on the 
formation of a two-electron dark state. We discuss the conditions under which such a state forms 
and describe the signatures that it leaves in transport properties such as the differential conductance 
and shotnoise. 
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I. INTRODUCTION 

Dark states (DSs) are a quantum-mechanical phe- 
nomenon originally discovered as a dark line in the flu- 
orescence of sodium atomsi. For a particular configura- 
tion of atomic transitions and driving fieldsii^>^, relax- 
ation can drive the atomic electron into a superposition 
state that is completely decoupled from the light field — 
this is the dark-state. 

The concept of the DS has been generaHzed to meso- 
scopic transporli^i^iii^i^ii^. In contrast to quantum 
optics, such systems are connected to electron reser- 
voirs, and enable one to study the infiuence of DSs 
on non-equilibrium transport properties such as current, 
shotnoiseii, and full counting statistical^. The first 
transport DS was proposed in a system similar to those of 
quantum optics, in which microwave fields were used to 
create a DS in a double quantum dol4>^. Subsequently, an 
all-electronic mechanism for the creation of DS has been 
described^i^ii^, in which coherent tunneling of electrons 
plays the role of the (classical) driving fields. The anal- 
ysed system was a triple quantum dot (TQD) in a trian- 
gular geometry in the strong Coulomb blockade regime. 

In the afore-mentioned systems, the DS is a single- 
electron state formed by interaction with the environ- 
ment. In the current work, we investigate the effects of a 
second electron on DS formation. This we do in a meso- 
scopic transport context, where a change of the chemical 
potential of the reservoirs can simply lead to the inclusion 
of two-electron states in the transport window. In par- 
ticular we study the infiuence of two-electron states on 
the transport through a TQD. The inclusion of doubly- 
charged states is particularly important because the for- 
mation of the original single-electron DS depends not 
just on destructive interference, but also on the strong 
Coulomb Blockade. The inclusion of two-electron states 
might therefore be expected to inhibit the appearance of 
DSs in such systems. However, as we will show, dark- 
state formation is in fact possible with two electrons, but 
only under certain circumstances. Even if the formation 
of the two-electron DS is incomplete, partial dark states 
still leave obvious signatures in the transport properties 
such as negative differential conductance and superPois- 
sonian shotnoise. 



This work comes against a backdrop of growing theo- 
retical interest in the transport properties of triple quan- 
tum dot a^^i^^i^^i^^ , and in the suppression of current 
due to interference phenomen a^^i^^i^^ as distinct from 
other current blockade mechanisms in quantum dots such 
as Coulomb^fi, spinal, isospin^^ and Franck-Condon^^ 
blockades. This work also has experimental relevance as 
a number of groups have published results of transport 
measurements on TQDs— i^^i^iZLiS^. The finite-bias cal- 
culations that we present here should facilitate the ex- 
perimental investigation of dark-state effects such as the 
break up of Coulomb blockade diamonds in the stability 
diagram of the TQD due to one- and two-electron DSs. 



II. MODEL 



Our system consists of three quantum dots (QDs) ar- 
ranged in a triangular geometry with a single relevant 
orbital level in each dot (see Fig.[l|). QDl and QD2 are 
connected to electron source reservoirs and QDS is con- 
nected to a drain reservoir. In the infinite bias Hmit, the 
rate at which electrons enter and leave the TQD is F, 
which we assume the same for all three leads. The levels 
in QDl and QD2 are coupled coherently to QDS with a 
tunnel amplitude Tc- 

We assume that the system is in the Coulomb block- 
ade regime, and we adjust the lead chemical potentials 
such that the only relevant charge states have zero, single 
and double excess electrons. We incorporate the single- 
electron charging energies into the energies of single elec- 
tron orbital levels. With the addition of a second elec- 
tron we associate an additional charging energy Uij, with 
i < j = 1,2,3, describing the locations of the two elec- 
trons. 
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(a) Two electrons in the same dot 




(b) Two electrons in different QDs 



Figure 1: The triple quantum dot (TQD) in triangular ge- 
ometry is connected to two sources and one drain. In the 
infinite bias limit, the rate at which electrons enter and leave 
the TQD is F, the same for all three leads. The QDs are 
coupled to each other by the coherent tunnel amplitudes Tc 
as shown. Zero, one or two electrons are allowed to be in the 
TQD at a time. The charging energy between two electrons 
in the same QD e.g. (a) U22 is, in general, larger than for two 
electrons in different QDs, e.g. (b) C/13. 



A. Hamiltonian 

The system Hamiltonian Hd of the closed TQD is a 
Hubbard-type model with three sites 

^^D =^E,n,„ + TcY^{d\^dy,„ + dlj:i„ + H.c.)+ 

i{7 (T 

3 3 

'^Uiihi^n^i+ ^ '^UijU^cfijc' , (1) 

where die is the annihilation operator, and n^o- is a cor- 
responding number operator of an electron in dot i with 
spin a. We assume spin-independent energy levels and 
denote the energy of the single-electron level of the QD 
i as Ei. Separating the Hamiltonian Hd into one-and 
two-electron parts, 

Hd = H^ + Hi + iJfi + Hii + H]i, (2) 

we explicitly construct the corresponding matrices in a 
basis of relevant many-body states. With the single- 
electron basis {|lcr), |2cr), |3cr)}, the Hamiltonians H^ for 
a single spin-up electron in the TQD and H^ for a spin- 



down electron in the TQD are 

[El Tc\ 
Hi=Hi=\ E2 Tc\. (3) 
\Tc Tc E^J 

The parallel-spin Hamiltonians H-^^ and H^ in the two- 
electron basis {|1ct2ct), |1ct3ct), |2ct3ct)} read 

(E1 + E2 + U12 Tc -Tc \ 

Hn=Hyi=\ Tc E1+E3 + U13 .(4) 

V -Tc E2 + E3 + U23J 

The opposite-spin two-electron Hamiltonian H-^i is a 9x9 
matrix, which is given in Appendix [Aj We denote the ba- 
sis of the above matrices the "localized basis". In subse- 
quent calculations we set all Ei = Eo + i?Gate = — eVoate 
(in the infinite bias calculations we set i?Gate = 0). 

A key parameter in our discussion of the two-particle 
dark state will be the difference between the charging 
energy Uu and the charging energy U12: 

Su = C/11 - C/12, (5) 

which we will call the "charging-energy difference". Since 
two electrons in a single QD are closer than two electrons 
in different QDs, one expects that the charging energies 
between electrons in different QDs {Uij; i ^ j) will be 
smaller than those for electrons in the same QD (Uu) 
and thus we focus here on Sjj > 0, (although note, for 
example, Refi^). We will discuss two sets of charging en- 
ergies: Firstly a highly symmetric situation where all Uu 
are equal and all Uij for i ^ j are also equal. Although 
this is the simplest situation, its high symmetry leads to 
non-generic features as we will see. We therefore consider 
a second set of charging energies with Un — U22 7^ C/33 
and U12 ^ C/13 = U23, which breaks the symmetry and 
leads to more typical results. 

The TQD is connected to three electron reservoirs that 
are described with the Hamiltonian 

^res = Y^^kc\k^Cika, (6) 

1=1,2,3 fc,cr 

where i labels the reservoirs (1,2 = source, 3 =drain). 
The TQD and the reservoirs are connected by the tunnel 
Hamiltonian 

i?T= E E^'^^.'^^ + H-c. (7) 

1=1,2,3 k,ij 

We assume spin-independent reservoir energies eik and 
tunneling amplitudes Vi^. 

B. Method 

We use two different approaches to calculate the cur- 
rent and the Fano factor of the TQD based on the 
sequential-tunneling (i.e. second order in tunneling 
Hamiltonian, Ht) Master equatio n^°i^^i^^ . The two ap- 
proaches are: 
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1. The generalized master equation 

In the infinite-bias limit, the second-order Born- 
Markov master equation is beheved to be exact for flat 
bands provided that coherences between system states 
are expHcitly included^i^. In the locaHsed basis, this ap- 
proach leads to the generalized master equation (GME). 
We assume that the chemical potentials are positioned 
far from other relevant energies, \fii\ > Uij,Tc,Vik,Ei, 
but such as to exclude three-electron states and states 
in which an electron occupies an excited orbital state of 
the QDs. We use the TQD Hamiltonian Hd and approxi- 
mate the chemical potential of the sources to be /.(i,2 = oo 
and the chemical potential of the drain to /is = — oo. In 
the GME approach, the density matrix p{t) contains en- 
tries for the populations of the empty state, and one- and 
two-electron states in the localized basis. Furthermore, 
p{t) also contains all coherences within each charge sec- 
tor. Within the GME, the time-evolution of the density 
matrix is then given by the Lindblad form^ 



'H^3'^ i^sapdl^-^dl^ds^p-^pdl^dsa) ■ (8) 



The second quantized operators appearing in this equa- 
tion are given in the relevant many-body basis in Ap- 
pendix [b1 We assume energy- and spin- independent tun- 
nel rates 



(9) 



with i = 1,2, 3. In the following we set all Fio- = F. The 
Lindblad master equation can then be written in super- 
operator formalism where the density matrix p is written 
as a vector. The equation of motion is then 



Pit) = Mp{t), 



(10) 



where M is the Lindblad super-operator. 

The stationary density matrix of the system is po = 
limt^oo p{t) and in practice this is determined as the 
nullvector of M, which is unique here. For our TQD, the 
steady-state current is then given by 

(/) J2 Ta. ((3a|po|3a) + 2(3? 3i|po|3t 3i) 

+ {la3a\po\la3a) + {2a3a\po\2a3a) 

+ (la3a|po|lcr3CT) + (2ct3ct|po|2ct3ct)) . (11) 



for the populations of the system eigenstates^S. We de- 
scribe this approach as the diagonalised master equation 
(DME). Unlike the GME, this approach is not restricted 
to the infinite bias limit, but the disadvantage is that 
some of the effects of coherence on the transport dynam- 
ics are lost. In the DME approach the TQD Hamiltonian 
Hd is first written in diagonal form 



N 



(12) 



fc=0 



where is the energy of eigenstate and = 21 

is the dimension of Hd without the empty state. The 
density matrix pnit) in the DME approach contains only 
the populations of states \Pk) and its time-evolution is 
given by the rate equation 

PD{t) = WpD{t), 

with elements of the rate matrix given by 



(13) 



Wk'k = Y.^^ 



[f{Ak^k+P^)m\dl\Pk') 
fI^)m\d^a\M\^ 



+ f{Ak'k 

N 

-^(/(A,,-/.,)i(Ai4i/3fc)r 

; 

where dia is the tunnel operator of the previous section, 
f{x) = (1 -I- e^f'^BT-^-i -g ^Yie Fermi function (with the 
temperature T and the Boltzmann constant fc^), and 
Ak'k = Efc' — Efc is the Bohr frequency of the transition 
from \(3k) to |/3^). 

In this case, the steady-state current is given by 

N 

(/)==eF3. [fi^k'k + P3)m\dl\M\HpDo)k 



k'k=0 



/(A,a-M3)|(/3feM3.|/3fc')l'(PDo)fc 



(14) 



with (pDo)fc as the fcth element of the diagonalized 
steady-state density matrix. 

In the infinite bias limit, the Fermi functions 

-Aij) = l 
J = 1,2; 

and 
In this 



of the sources leads are lim^^^oo /(Afc/fc 



lim 



/^j — *oo 



/(Afe.fc+Mj) = 0, 



and 

and lim^3^_oo /(Afe'fc -I- /i3) = l, 

lim^3^_oo /(Afe'fc - /i3)=0 for the drain, 
limit, the current is then 



N 

■ E 

k,k'=0 



m\di\M\\pD,)k. (15) 



S. Rate equation in the diagonalised basis 



III. STATIONARY TRANSPORT IN THE 
INFINITE BIAS LIMIT 



An alternative approach is to first diagonalise the sys- 
tem Hamiltonian and then write down a rate equation 



We first consider the infinite bias limit. It is in this 
limit that previous calculations on the single-electron 
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Figure 2: Stationary current (/) through the TQD as a func- 
tion of the charging-energy difference scaled with the tun- 
nel amplitude Su /Tc for different ratios of Tc/T and for two 
choices of charging energies. Results are shown for both GME 
and DME calculations (the DME calculations are independent 
of the ratio Tc/T) and both calculations show that at Su = 0, 
the current is zero. This is attributed to the formation of the 



two-electron dark state. In Fig. 2(a) (symmetric charging en- 
ergies) the GME calculation shows an antiresonance of finite 
width about the DS, whereas in the DME calculation the DS 
appears as a point of discontinuity. For asymmetric charging 
energies, Fig. 2(b) both calculations show an antiresonance 
of finite width about the DS. 



TQD have been performed^i^ii^ and where we expect the 
GME to be exact. We therefore discuss the GME results 
first and then return to a comparison of GME and DME 
calculations in this Hmit. 
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Figure 3: Gap width in Su/Tc of the antiresonance at the DS 
as a function of tunnel rate per tunnel ampli tude V/Tg for 
the GME approach with the parameters of Fig. 2(a) The size 
of the gap is given by the width of the antiresonance when 
the current has half the size of the maximum ^^"^^f^^ . 



charging-energy difference vanishes {6u = 0) the current 
is zero. At this point the system is trapped in the two- 
electron DS 



Pdark 



with 



5« 



(16) 



(17) 



For Su = 0, I'Z'o) is an exact eigenstate of the two-electron 
Hamiltonian. As since this state has no occupation on 
QD3 the electrons in the TQD cannot leave the TQD to 
the drain and no further electrons can enter the TQD 
due to the Coulomb blockade, this state is completely 
dark with zero current. This situation is analogous to 
the single-electron DS but here we have two electrons 
in forming a product state of a spin-up and spin-down 
single-particle DSs. 

The current in the GME approach for Su has a 
dark resonance profile with a clear antiresonance. The 
width of this antiresonance as a function of the tunnel 
rate is shown in Fig. [31 For small ratios of T/Tc < 1 the 
gap increases Hnearly from T/Tc = (corresponding to 
the limit of Tc 00 where the DME current and the 
GME current coincide). For higher ratios of T/Tc the 
size of the gap reaches a maximum at T/Tc « 3.7, and 
then decreases. 



A. Stationary Current 

Figure E] shows the stationary current (/) as a func- 
tion of the charging-energy difference Su , which we scale 
with the tunnel amplitude Tc. At the point where the 



B. Shotnoise 

In the infinite-bias limit, we can write down the n- 
resolved master equation 



(18) 
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(a) U22 = t^33 = Uu and I/13 = U23 = U12 




(b) C/11 = U22, U12/TC = 10, U33/TC = 15 and 
Uw/Tc = U23/TC = 11 

Figure 4: The Fano factor F{0) as a function of the charging- 
energy difference with results for both GME and DME cal- 
culations. Fr om the GME calculati on for both symmet- 
ric (Fig. |4(a) I and asymmetric (Fig. 4(b) I charging ener- 
gies, we observe that the shotnoise is highly super-Poissonian 
(F(0) > 1) both near the dark state Su/Tc ~ and in 
the limit of Su/Tc ^ 1. In the DME approach, only the 
asymmetric case shows super-Poissoni an be haviour near the 
dark state. The Fano factor in Fig. 4(a) is Poissonian for 
Su/Tc < 1. 



in the infinite-bias limit of the DME. 

The advantage of the n-resolved master equation is 
that it can be used calculate not only the stationary cur- 
rent (given by (/) = e Tr(Mjpo) in this language), but 
also the full counting statistics of the currenli^. Here 
we concentrate on the shotnoise^i, and from the method 
used in Ref.— , the zero-frequency Fano factor i^(0) can 
be expressed as 



(21) 



k=l 



with Ck = Ea=0 T,'^",,,=oiMj)af3Vpk{vkf) ^ {M j) ^ ^,{pq) ^ 

and \k the eigenenergies (Ao = 0). Here, {Mj)ai3 are the 
components of the jump operator, vpk are components 
of the matrix of eigenvectors of the Liouvillian M, and 
(po)^! is the ^th component of the steady state density 
matrix. is the number of populations in the density 
matrix and Nj^i is the dimension of M . In our GME 
calculation was A^a/ = 117 and in the DME Nm = 21 
and in both cases = 21. 

FigurelHshows the Fano factor F(0) as a function of the 
charging-energy diflFerence Sy /Tq. Two zones of super- 
Poissonian behaviour are observed: in the region near 
the dark state <5;7 = and also for large charging-energy 
diflerence &u /Tq ^ 1 . Enhancement of the Fano factor 
near the dark state is expected and can be understood 
as a dynamical channel blockade^ with the dark state 
playing the role of the weakly coupled channel. The value 
of the Fano factor near J(7 = is w 6, which is roughly 
twice that found for the single-electron dark stated. The 
precise mechanism behind the highly super-Poissonian 
behaviour at large 5u /Tc is not yet clear. However, it 
is a quantum coherent phenomenon, as the Fano factor 
drops markedly when dephasing is included^. 



C. Difference between GME and DME approach 



for p'-"-' , the partial density matrix of the system after n 
electrons have passed though the TQD^. Here the total 
LiouvilHan M has been decomposed into two parts: A/q, 
which describes the evolution of system without electron 
transfer to the drain, and Mj, the "jump operator", which 
transfers an electron to the drain. The matrix elements 
of the jump operator are obtained from 



3<J 



=U 



in the GME approach, and as 



(19) 



(20) 



Fo r asy mmetric charging energ ies, the current, 
(Fig. |2(b)| and the Fano factor (Fig. |4(b)[ ) of the DME 
calculation show the same general features as the GME 
results, with qualitative differences becoming more evi- 
dent for s mall values of Tc/T. In the symmetric case 
(Fig. 2(a) ), however, a qualitative difference between the 
GME and DME results is observed in the region of the 
DS. Both calculations show a complete current suppres- 
sion at i5;7 = but, whereas the GME shows a broad 
antiresonance about the point Sjj = 0, the DME shows a 
discontinuity at this point. The discontinuous behaviour 
of the symmetric DME current is a reflection of the in- 
creased degeneracy of the TQD eigenstates when all Uij 
are equal. At this point, the dark state resides in a 
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degenerate subspace with two other orthogonal states 
1 



1 



lT3i)-|2t3i) + |3TU) + |3T2i)) (22) 
-|iTU)-|iT2i)-|2T2i)-|2Tii» 

(23) 



+ ^|3T3i). 

Due to this degeneracy, one has, in principle, a free- 
dom of choice as to which combination of these states 
to treat as the eigenstates to be entered into the DME 
machinery. Different choices result in different results. 



In obtaining the plots of Figs. 2(a) and 4(a) we chose 
the combinations as above since, from the GME analy- 
sis, we expect a DS to form. However, without this a 
priori knowledge, any linear combination of the above 
vectors appears as good as any other, and if any combi- 
nation other than that given explicitly above is chosen, 
the dark state will not be observed since all three eigen- 
states will have a finite population on QD3. We believe 
that this undesirable feature of the DME should be re- 
moved by a more complete second-order master equation 
treatment^'^i^'^i^, in which the coherences are han- 
dled properly. For Uu ^ Ujj, (i ^ j) or for asymmetric 
charging energies as in Fig. 2(b) this degeneracy is lifted, 
the choice of eigenstates is unique and the above prob- 
lems do not occur. 

The Fano factor for symmetric charging energies cal- 
culated with the DME in Fig. |4(a) is also interesting be- 
cause, although it reproduces the high values at Su/Tc ^ 
1, its behaviour is Poissonian, F(0) = 1, at small values 
of 5ij /Tc where the GME shows super-Poissonian be- 
haviour. It appears then that, for symmetric charging 
energies the effects of the two-electron dark state in the 
in the DME approach are restricted solely to the singular 
point 5u/Tc = 0. 



IV. FINITE BIAS 

As long as the temperature is low and all energy levels 
of the TQD are well within the transport window, the re- 
sults of the GME in the infinite bias limit are reliable^^*^. 
However, the GME is not applicable away from this sit- 
uation, and in this case we employ solely the DME. Our 
finite-bias calculations allow us to construct stability di- 
agrams for the system and determine the implications 
of the one- and two-particle dark states for tunnel spec- 
troscopy measurements^. 

We first consider the situation with 6ij — ^ such that 
the dark state is an eigenstate of the TQD Hamilto- 
nian. In accordance with the above discussion, we choose 
asymmetric charging energies such that methodical prob- 
lems are not an issue. The resulting current is shown in 
Fig. [5] and the differential conductance in Fig. [H The 
bias voltage Veias is chosen symmetric to Voate (hence 
M; = and /i; = — VBias/2). For negative bias volt- 

age, Veias < 0, the electrons enter the system through 




-0.6-0.4-0.2 0.2 0.4 0.6 O.S 1 1.2 1.4 
VG.te(mV) 

Figure 5: Stability diagram for the TQD with parameters 
such that both one- and two- electron dark states can form; 
in particular, Su — 0. Here the current is plotted as a function 
of the source-drain bias voltage Vsias and gate voltage Voate- 
For positive bias, most electrons enter the system on the left 
side where two QDs couple to the environment, and leave to 
the right side where one QD is coupled to the environment. 
Only for this bias direction can the DSs occurs. The DSs 
are visible as breaks in the borders of the Coulomb diamonds 
(marked as DSl for the one electron dark state and as DS2 
for the two electron dark state); the gap near the center is 
due to the single-particle dark-state, that to the right, the 
two-particle dark state. The parameters used are: Uu = 
U22 = U12 = ImeV, U33 = 1.2meV, U13 = U23 = 0.95meV, 
tunnel amplitude Tc = O.lmeV, tunnel rate F = lO^eV and 
temperature T = 150mK. 



QDS and leave it through QDl or QD2 and the lower 
halves of Figs. [5] and [6] show the familiar diamonds of the 
Coulomb blockade. The light-colored lines in the differ- 
ential conductance diagram that lie outside the diamonds 
correspond to excitation energies of the TQD. 

For VBias > 0, this picture changes drastically and the 
stability diagram is dominated by the effects of the DSs. 
Both one- and two-electron DSs are visible. In the cur- 
rent diagram they show up as breaks in the Coulomb dia- 
monds, and in the differential conductance they give rise 
to clear lines of negative differential conductance (NDC). 
Note that, even for large Veias > 0, it is possible that no 
current fiows. This strong asymmetry between Veias < 
and Veias > demonstrates that the dark states can lead 
to a strong rectification of the current. For small charg- 
ing detuning Su ^ 0, the stability- and the conductance- 
diagrams look similar to Fig. [5] and Fig. [6] but with a 
small current fiowing along the line DS2. 

It is also of interest to consider what happens in 
this system far from the dark resonance. A differential 
conductance diagram for this situation with symmetric 
charging energies is shown in Fig. \7\ Note that in this 
figure, we are far from the DS such that the DME cal- 
culation should be reliable, even for symmetric charging 
energies. With this set of parameters, the one-electron 
DS can still form and this is clearly observed as a break 
in the Coulomb diamond and a diagonal line of nega- 
tive differential conductance. Moreover, traces of a two- 
electron dark states are still visible, manifest as lines of 
negative differential conductance in the two-electron sec- 
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This state is similar to the dark state of Eq. fT7|) but 
with only contributions with two electrons in each dot. 
With an increasing difference between Ua and Uij for 
j ^ i the line of negative differential conductance which 
corresponds to NC2 becomes stronger and also the weight 
of the state Eq. (|24|) increases. 



CONCLUSIONS 



Figure 6: Differential conductance as a function of bias volt- 
age Veias and gate voltage Voate with the same parameters as 
Fig.O For Veias > 0, two lines of negative differential conduc- 
tance (NC) are visible, which can be attributed to the one- 
and two-electron DS. When a DS is hit the current decreases 
to zero and the differential conductance becomes negative. As 
in the current, the diamond pattern is broken by the forma- 
tion of DSs. 
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Figure 7: As Fig.[6l but with finite charging-energy difference 
&JJ = 0.5meV and symmetric charging energies. Despite the 
fact that the genuine two-electron dark state (Eq. I|17p l can 
no longer form, lines of negative differential conductance are 
seen in positions similar to that in Fig. [6l A strong new line 
of NDC is also observed (marked NC2) which is related to the 
state of Eq. I|24p . More structure is visible here as compared 
with the &u plot since degeneracies have been lifted. The 
parameters used here were Un = ImeV, Uij — 0.5meV for 
j ^i,Tc = lOOAteV, r = 10/xeV and T = lOOmK. 



tor. These occur because, although the full dark-state is 
no longer an eigenstate of the system Hamiltonian, there 
still exists related states which, although not completely 
dark, can only support a comparatively small current. 
In particular, there is one NDC line (marked NC2) in 
the two-electron region of the stability diagram for which 
the suppression is particularly strong. The corresponding 
stationary density matrix has very low occupation on the 
third QD and predominantly given by pstat ~ Wd2){'I^D2\ 
with 



D21 



/2 V IT U 



V2 



4^4 J io). 



(24) 



We have shown in this article that the TQD in trans- 
port can exhibit current suppression due to a two-particle 
DS. This represents an extension of the single-electron 
DS concept famihar from optics to the multiple-electron 
domain. In the TQD, the two-electron DS can ex- 
ist only under the condition that the charging energies 
Un ~ U22 = U12 are equal. Importantly though, even if 
these energies are not equal, traces of the DS still remain, 
and should be observable in experiment. 

Specifically, these traces are a large Fano factor near 
the DS and pronounced lines of negative differential 
conductance in the stability diagrams. We note that an 
optical two-electron DS was presented in Ref4i. An anal- 
ogous mesoscopic transport configuration would not lead 
to a DS formation and hence current blockade, since the 
dark-state has finite occupation of all states in the local- 
ized basis. 

On a technical level, we have compared the GME and 
DME approach in the infinite bias Hmit. We believe that 
the gap between these two methods can be effectively 
bridged with a more complete second-order master equa- 
tion treatment, e.g.— >^>^>^>^. This should enable us to 
calculate transport properties at finite bias without hav- 
ing to concern ourselves with the discontinuity exhibited 
by the DME. 

Further work includes investigation of dark states with 
still higher electron numbers and the infiuence of mag- 
netic field. 
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Appendix A: OPPOSITE-SPIN TWO-ELECTRON 
BASIS 



In the basis { | IT li) , | IT 2i) , | IT 31) , 1 2T 2i) , 1 2T li) , 1 2T 3i) , 
|3T3|), |3Tli), |3T2T)}, the Hamiltonian for two electrons 
of opposite spins reads 



8 



i+Uii 





Tr 














Tr 








-C/1 ~h-C/9~hcy 1 9 

-•—'X 1 ' ^ 1 ^ X 


Tc 

















Tc 


Tc 


Tc 













Tc 

















2E2+U22 





Tc 








Tc 














E2+EI+U12 


Tc 





Tc 














Tc 


Tc 


E2+E3+U23 


Tc 














Tc 








Tc 


2E3+U33 


Tc 


Tc 


Tc 











Tc 





Tc 


E3+E1+U 13 








Tc 





Tc 








Tc 





E3+E2-' 



(Al) 



Appendix B: TUNNEL OPERATORS 

The second-quantized operators appearing in Eq.® can be written in the relevant many-body basis as: 



rfsT - |0)(3t| - |1T)(1T3T| - |2t)(2T3T| + |U)(3TU| + |2i)(3T2i| + |3i)(3T3i|, (Bl) 

d3l = |0)(3i| - |U)(lT3i| - |2t)(2t3i| - |3t)(3t3i| - |U)(U3i| - |2i)(2i3i|, (B2) 

= |lt)(0| + |iT2T)(2t| + |1T3T)(3T| + liTUXUI + |iT2i)(2i| + |lT3i)(3i| (B3) 

4^ = |2T)(0| - |1T2T>(1T| + |2T3T)(3T| + |2TU)(U| + |2T2i)(2i| + |2T3i)(3i|, (B4) 

= |U)(0| - liTUXiTI - |2TU>(2TI - |3TU)(3t| + |U2i)(2i| + |U3i)(3i|, (B5) 

4^ = |2i)(0| - |lT2i)(lT| - |2T2i)(2T| - |3T2i)(3T| - |U2i)(U| + |2i3i)(3i|. (B6) 
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